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Abstract 

We measure the statistical polarization of quadrupolar nuclear spins in a sub-micrometer 
(0.6 /xm^) particle of GaAs using magnetic resonance force microscopy. The crystalline sample 
is cut out of a GaAs wafer and attached to a micro-mechanical cantilever force sensor using a 
focused ion beam technique. Nuclear magnetic resonance is demonstrated on ensembles contain- 
ing less than 5 x 10^ nuclear spins and occupying a volume of around (300 nm)'^ in GaAs with 
reduced volumes possible in future experiments. We discuss how the further reduction of this de- 
tection volume will bring the spin ensemble into a regime where random spin fluctuations, rather 
than Boltzmann polarization, dominate its dynamics. The detection of statistical polarization in 
GaAs therefore represents an important first step toward 3D magnetic resonance imaging of III-V 
materials on the nanometer-scale. 

PACS numbers: 76.70.-r, 05.40.-a, 68.37.Rt, 85.85-|-j 



I. INTRODUCTION 



Recent years have seen the development of a wide range of semi-conducting nanostruc- 
tures including quantum wells (QWs), nanowires (NWs), and quantum dots (QDs). Re- 
searchers have devoted particular attention to making devices from III-V materials such as 
GaAs, whose high electron mobility and direct band gap make it a critical component of 
today's semiconductor technology. 111-V systems are extremely versatile in large part due to 
techniques such as molecular beam epitaxy (MBE) and metal-organic chemical vapor depo- 
sition (MOCVD), which enable the growth of complex heterostructures with nearly perfect 
crystalline interfaces. As a result, applications range from integrated circuits operating at 
microwave frequencies to light-emitting and laser diodes to quantum structures used for 
basic research. 

While a variety of techniques exist to characterize and image these nanostructures, in- 
cluding scanning electron microscopy (SEM), tunneling electron microscopy (TEM), and 
x-ray crystallography, so far it has been impossible to measure single nano-structures using 
magnetic resonance imaging (MRI). In larger structures, MRI is a powerful technique allow- 
ing for the three-dimensional (3D), sub-surface imaging of the density of particular nuclear 
magnetic moments. However, conventional nuclear magnetic resonance (NMR) techniques, 
in which the spin signal is detected by an inductive pick-up coil, are limited to detection 
volumes of several [xm on a side or larger.^ A net polarization of at least 10^^ nuclear spins 
is typically needed to generate a detectable signal; nanometer-scale samples simply do not 
contain enough spins to be detected. In the past few years, a more sensitive force-detected 
version of MRI has been demonstrated on nanometer-scale samples.^ Using magnetic res- 
onance force microscopy (MRFM) to measure the statistical polarization of spin-1/2 ^H, 
Degen et al. made 3D images of single virus particles with a resolution better than 10 nm.^ 

Here we take a step toward applying this technique to quadrupolar (spin-3/2) nuclei, 
specifically Ga and As, in a nanometer-scale particle. We demonstrate the detection of 
statistical polarizations of ^^Ga, ''^Ga, and ^^As in a 0.6 |J.m^ particle of crystalline GaAs. 
The mechanical detection of NMR in GaAs was first demonstrated in 2002 by Verhagen et al. 
and Thurber et al.^'^ The smallest reported detection volume of 600 |J.m'^ contained more than 
10^^ nuclear moments.^ In 2004, Garner et al. reported forced-detected NMR signal from 
10^° moments in a GaAs wafer. Those experiments measured either the thermal equilibrium 
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FIG. 1: SEM micrograph of the Si cantilever with a GaAs sample attached. (A) shows the cantilever 
protruding from a Si chip. (B) clarifies the geometry showing the paddle and mass-loaded end of 
the cantilever, and (C) is a detailed view of the tip of the mass-loaded cantilever with the GaAs 
sample attached. A layer of Pt is visible at the very tip of the GaAs particle. 

polarization or an optically enhanced polarization of Ga and As spins. Our experiment has 
a detection volume of about 0.03 (xm^ pa (300 nm)^ equivalent to less than 5 x 10^ spins of 
any one of the constituent isotopes. Such a volume is far too tiny to detect via conventional, 
inductively-detected magnetic resonance, although the number of spins is not yet small 
enough that its polarization is dominated by statistical fluctuations. Future reductions in 
detection volume, however, will enter this regime and will require techniques like the one 
demonstrated here. 

II. BOLTZMANN VS. STATISTICAL POLARIZATION 

Conventional magnetic resonance signals originate from the mean polarization of nuclear 
spins in an external magnetic field - the so-called Boltzmann polarization. Although this 
polarization is quite small, it dominates the spin signal for large ensembles of nuclear spins. 
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As the size of the spin ensemble decreases, the amphtude of the polarization fluctuations 
eventually exceeds the mean polarization.^ This variance, sometimes called the statistical 
polarization, then becomes a more useful signal for MRI than the mean polarization. 

Statistical polarization arises from the incomplete cancellation of randomly oriented spins. 
For any given direction, the net polarization can be either positive or negative and will fluc- 
tuate on a time scale that depends on the flip rate of the spins. Several MRFM experiments 
have detected statistical polarizations of spin-1/2 nuclear spins^"^^ and demonstrated their 
use for nanometer-scale nuclear MRI.'^'^^ 

In order to understand the regimes in which either Boltzmann or statistical polarization is 
important, consider an ensemble of spins with spin quantum number J. The Hamiltonian 
of a single spin in the presence of a magnetic field B along the z is H = —fizB = —h'-yBIz 
where fiz is the magnetic dipole moment operator along z, h is Planck's constant, 7 is the 
gyromagnetic ratio, and is the nuclear spin angular momentum operator along z. Sta- 
tistical mechanics predicts the equilibrium distribution at a temperature T to be a Boltz- 

_ H 

mann distribution. The partition function Z = Tr{e ''bt^ contains all the information 

ft 

about the nuclear spin polarization in the system. The density matrix p = -^e of 
the spin system can be used to compute the mean = NTr{fizp} and the variance 
a\j^ = N{Tr{filp} — (Tr{/iz/)})^) of the ensemble's magnetization along z. Even at cryo- 
genic temperatures (T ~ 1 K) and high magnetic fields {B ~ 10 T), h'-fB ^ ksT. Therefore, 
keeping terms only up to first order in we find: 

^k-N'^ih-rf. (2) 

Since Nh'yl corresponds to 100% spin polarization, one can define the statistical nuclear 
polarization as SNP = = \J ^jfjf and the Boltzmann nuclear polarization as BNP = 
-^f^ = Statistical polarization dominates the system {SNP > BNP) when the 

number of spins in an ensemble is less than the critical number, 

Equivalently, for a material with a nuclear spin density na, where n is the number density of 
the nuclear element and a is the natural abundance of the isotope of interest, one can define 



4 



a critical volume K = For volumes smaller than V^, magnetic resonance experiments 
should be designed to detect SNP rather than BNP. 

III. MRFM TECHNIQUE AND APPARATUS 

We measure the presence of a particular nuclear isotope using an MRFM protocol which 
cyclically inverts statistical spin polarization.^^ In a magnetic field Btotah the frequency 
of a transverse RF magnetic field Bi is swept through the nuclear resonance condition, 
/rf = ^-Btotai- If done adiabatically, this sweep induces the nuclear spins to invert - a 
process known as adiabatic rapid passage. In the strong spatial magnetic field gradient 
near a magnetic tip, these inversions produce a time- dependent force. This force is in turn 
detected as the displacement of an ultrasensitive cantilever. 

Our MRFM experiment is carried out in a sample-on-cantilever configuration in which 
the sample is affixed to the end of a single-crystal Si cantilever^'' as shown in Fig. 1. We 
arrange the cantilever in a "pendulum" geometry such that the sample is positioned above 
a FeCo magnetic tip, depicted in Fig. 2. The magnetic tip, which is mounted on a separate 
chip, is patterned on top of an Au microwire, which acts as an RF magnetic field source and 
is shown in Fig. 3.^° 

The cantilever measures 120 |J.m x 4 |j,m x 0.1 |a.m and - loaded with the GaAs sample - 
has a mechanical resonance frequency /c = oJc/ (27r) = 3.7 kHz and an intrinsic quality factor 
Q = 4.0 X 10'^ at T = 1 K. By measuring the cantilever's thermal motion, we determine 
its effective spring constant to be A; = 120 |j,N/m. The MRFM apparatus is isolated from 
vibrational noise and is mounted in a vacuum chamber with a pressure below 10~^ mbar 
at the bottom of a ^He cryostat. The motion of the lever is detected using 100 nW of 
1550 nm laser light focused onto a 12 |^m-wide paddle and reflected back into an optical 
flber interferometer. The microwire used to produce the transverse RF magnetic fleld is 
2.5 [xm-long, 1 |j.m-wide, and 0.2 (xm-thick. The FeCo tip is shaped like a bar, sits on top 
of the microwire, and produces a spatially dependent field Btip(r). It has a top width of 
270 nm, a bottom width of 510 nm, a length of 1.2 [im, and a height of 265 nm as shown 
in Fig. 3. To make sure that the FeCo tip is fully magnetized along the z, an external 
magnetic field B = Bz is applied with B = 2.65 T. During measurements, the distance 
between the FeCo tip and the closest point on the sample is typically 100 nm such that the 
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FIG. 2: Representation of the MRFM apparatus at the bottom of the cryostat. The microwire is 
shown in yellow, the FeCo tip in blue, and the GaAs sample in red, with a Pt layer on its end in 
white. The green region above the FeCo tip depicts the resonant slice during measurement. B, 
the cantilever shaft, and the magnetization of the FeCo tip are aligned along z. Near the FeCo 
tip, current flows in the wire along y, while the lever displacement and Bi at the position of the 
sample are directed along x. 

static magnetic field gradient relevant to MRFM is on the order of 5 x 10^ T/m, where 

Btotai = B + Btip and x is the direction of cantilever oscillation. The maximum |Btip| for 
this spacing is about 0.1 T. Smaller spacings can result is larger and |Btip|, although 

they also result in larger measurement noise. Interactions between the magnetic tip and the 
sample at such small gaps, known as non-contact friction, lead to mechanical dissipation 
in the cantilever. ^^'^^ In our experiments, these effects reduce the quality factor Q of the 
cantilever to 1.0 x 10^. In addition, we damp Q down to ~ 400 using active electronic 
feedback. Given the narrow natural bandwidth of our high-Q cantilever, damping serves 
to increase the bandwidth of our force detection with out sacrificing signal-to-noise ratio 
(SNR).2i 
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FIG. 3: SEM micrograph of the Au microwire with integrated FeCo tip. The structure is patterned 
on a Si chip. 



If the rate of the RF frequency sweeps used to invert the nuclear spins is slow enough 
and the amplitude of Bi large enough, the initial population distribution among the nuclear 
spin energy levels is completely inverted. As a result, the net magnetization is made to 
flip along the z. For spin-1/2 nuclei, the criterion for adiabatic inversion is given by a = 
2Tr'y^ Bf / (uc^yif) ^ 1, where f2RF/(2vr) is the amplitude of the frequency modulation around 
the center frequency /rf of the transverse RF field Bi}^ The criterion for quadrupolar 
nuclei, in general, is more complex. However, a complete inversion of the initial population 
distribution over the quadrupolar energy levels can be achieved in the limit of both a 3> 1 
and (3 = 'jBi/Qq ^ 1, where fiQ/(27r) is the quadrupolar frequency. The frequency sweeps 
used here are designed to meet these conditions and follow the form used in Poggio et al.^° 
Therefore, by sweeping through f^p at a frequency 2/c, we are able to modulate the spin 
polarization at fc- The resulting spin inversions produce a force that drives the cantilever 
at its resonance frequency. An ensemble of spins at position -r with a statistical variance in 
its 2;-magnetization a^^ produces a force on the cantilever with variance. 



Using our knowledge of the spring constant k, we determine ap by measuring the variance 
of the cantilever's oscillations on resonance. The correlation time of is limited by the 
relaxation rate of the spins in the rotating frame. 




(4) 
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FIG. 4: MRFM signal from the statistical polarization of ^H, ^^Ga, '^^Ga, and ^^As. Black dots 
show the resonant force variance ap as a function of the center frequency /rf- Solid red lines 
represent adjacent-averaging of the data as a guide to the eye. Inset is a zoomed-in view of the 
spin-3/2 nuclear resonances: ^^Ga, "^^Ga, and "^^As. Error bars represent the standard error of ap 
calculated as in Degen et al.^^ Data points represent 1400 s of averaging for ^^As, 600 s for ^^Ga, 
600 s for ■^iGa, and 300 s for ^H. 

IV. RECEPTIVITY IN MRFM OF STATISTICALLY POLARIZED ENSEMBLES 

Due to differences in magnetic moment and statistical polarization, two ensembles con- 
taining the same number of nuclei but each of a different isotope, produce different magneti- 
zation variances. This difference is contained within the concept of receptivity. Receptivity 
is a value defined for the purpose of comparing the expected NMR signal magnitudes for 
equal numbers of different nuclear isotopes. For MRFM of statistically polarized ensembles, 
we define a receptivity, -Rn.mrfm oc 7^/(1 + 1), where -Rn.mrfm is normalized to 1 for -"^H. The 
factor 7^/(/ + 1) is proportional to the magnetization variance expected from an ensemble 
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-Rn.mrfm 
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0.293 


0.466 


0.147 
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0.071 


0.142 


0.025 




-Rv.MRFM 




1 


0.056 


0.058 


0.046 








1 


0.013 


0.018 


0.008 



TABLE I: Properties relevant for a statistically polarized MRFM measurement. 



of spins as defined in (2). From (4) the variance, multiplied by the square of the magnetic 
field gradient, results in the resonant force variance measured in MRFM. 

On the other hand, conventional NMR collects an inductive signal due to a Boltzmann 
polarization. In this case, receptivity can be defined as -Rn.cohv oc 7^/(J + l), where -Rn.cohv is 
similarly normalized to 1 for ^H. Here the factor 72/(1+ 1) is proportional to the Boltzmann 
polarization as defined in equation (1). The remaining factor of 7 results from the fact 
that conventional NMR measures the inductive response of a pick-up coil to magnetization 
precessing at a frequency proportional to 7.2^ 

As can be noted in Table I, MRFM receptivity scales more favorably than conventional 
receptivity for low-7 nuclei such as those found in GaAs. In real experiments, comparisons 
are often made between signals from two different isotopes contained in the same volume. 
In the comparison of volumes rather than number of nuclei, one must also take into account 
the number density n of each element in the material and its natural isotopic abundance a. 
Volume receptivity therefore also includes the factors of n and a: -Rv,mrfm oc nai?N,MRFM 
and -Rv.conv oc nai?N,conv 
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V. MRFM MEASUREMENTS 



Here we study a sub-micron sized particle cut from the surface of a GaAs wafer. The 
GaAs sample is affixed to the cantilever tip using a focused ion beam (FIB) technique. First, 
a thin layer of Pt is deposited over a small area of a GaAs wafer to protect the sample from 
potential ion damage. Then, a lamella measuring 3 (xm x 2 [im x 0.3 |J.m is cut out from this 
area of the wafer. Next, the lamella is welded with Pt to a nearby micro-manipulator and 
transferred to the tip of an ultrasensitive Si cantilever. Finally, the particle is Pt-welded to 
the cantilever tip and cut to its final dimensions: 2.4 |am x 0.8 (xm x 0.3 [im = 0.6 \xm^ 
(Fig. 1). The side of the sample which formerly was part of the wafer surface is oriented 
such that it faces away from the cantilever. A roughly 200 nm-thick layer of the original Pt 
protection layer remains on this surface of the particle. Throughout this process, only the 
mass-loaded end of the cantilever is exposed to either the ion or electron beams. Special 
care is taken never to expose the cantilever shaft to in order to avoid structural damage 
or deposition of material on its surface. Even short exposure can lead to the permanent 
bending of the cantilever and a reduction of its mechanical Q. 

MRFM signal measured from this GaAs particle at -B = 2.65 T and temperature T = 1 
K is plotted as a function of the RF center frequency in Fig. 4. Resonances from all three 
isotopes (^^Ga, ^^Ga, and ^^As) in GaAs are visible. In addition, a strong resonance 
appears in the spectrum due to the thin layer of adsorbed hydrocarbons and water that 
coats surfaces which have been exposed to ordinary laboratory air.'^'^^ Each resonance is 
measured with the GaAs particle positioned at slightly different x and y positions in the 
vicinity of the FeCo tip. In each case, however, the spacing along z between the end of the 
particle and the top of the FeCo tip is 100 nm. Similar magnitudes of Bi are used in each 
case, which we quantify in the discussion of Fig. 5. The frequency modulation amplitude 
fiRF/(27r) is 400 kHz for ^H, 100 kHz for ^SGa and ^^Ga, and 50 kHz for ^^As. Each data 
point represents 300 s of averaging for ^H, 600 s for ^^Ga and ^^Ga, and 1400 s for ^^As. 
While the SNR for some peaks is small - ''^As and ^^Ga in particular - each peak is confirmed 
by at least one other experiment performed at a different magnetic field B. The appropriate 
shift in carrier frequency is observed in each case. 

The rotating-frame spin correlation time observed for ^H is on the order of 100 ms, 
which is consistent with spin correlation times measured in similar experiments on surface 
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hydrocarbon layers. ^'^^ For the quadrupolar isotopes in GaAs, we measure Tm to be around 
500 ms. Thurber et al. report on the order of several seconds in MRFM measurements 
of Boltzmann polarized quadrupolar spins in a GaAs wafer. ^ Given that the adiabaticity 
parameter a is similar to that used in our experiments, the difference in Tm is likely due to 
the large difference in magnetic field gradient in the two cases. Gradients in our experiments 
exceed 10^ T/m while gradients used by Thurber et al. are 100 times smaller. High magnetic 
field gradients on this order have been shown to limit in similar MRFM experiments.^^ In 
particular, mechanical noise originating from the thermal motion of the cantilever couples 
through strong magnetic field gradients to produce nuclear spin relaxation in the statistically 
polarized ensemble. 

The central frequency, amplitude and, width of the resonance peaks depend on the various 
experimental parameters including 7, B, the spatial dependence of i?tip, the shape of the 
sample, its position relative to the FeCo tip, and the form of the adiabatic sweep waveform. 
Roughly, however, one can say that the low-frequency onset of each resonance should occur 
at ^B. Note that the peak magnitudes in Fig. 4 do not scale with -Rv.mrfm since both the 
volume of the material detected and the magnetic field gradient in that volume are different 
for each measured peak. The striking difference in signal amplitude between the hydrocarbon 
layer and the quadrupolar nuclei in the GaAs particle is mostly due to the smaller gradients 
present inside the GaAs particle compared to those present at the hydrocarbon layer. Due to 
the 200 nm Pt layer covering the tip of the GaAs particle, the containing layer is 200 nm 
closer to the FeCo tip than any of the isotopes in GaAs. As a result, for the same sample tip- 
sample spacing, the nuclei experience about ten times higher than the quadrupolar 
isotopes - resulting in force variances 100 times larger from the same magnetization. We 
discuss the effect of each parameter on the resonances in more detail in Section VII. 

VI. NUTATION MEASUREMENTS 

Using the method described in Poggio et al.,^° we also measure the rotating-frame am- 
plitude of Bi. Pulses of variable length are inserted in the adiabatic sweep waveform every 
500 cantilever cycles (135 ms). The measured force variance in spin nutation experiments 
of ^^Ga is plotted in Fig. 5 for a spacing along z between sample and FeCo tip of 100 nm. 
The amplitude of the rotating RF magnetic field is inferred by fitting the data to a decay- 
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ing sinusoid. The frequency of the sinusoid corresponds to the Rabi frequency 7i?i of the 
isotope in question and the decay rate is related to the spatial inhomogeneity of i?tip within 
the detection volume. The measured Rabi frequency of 208 kHz for ''^Ga corresponds to 
Bi = 16 mT. This measurement represents the rotating-frame amplitude Bi in the region 
of the GaAs particle closest to the FeCo tip, where the gradients and the resulting contri- 
bution to the MRFM signal are largest. Similar nutation measurements done using the 
containing layer, which is 200 nm closer to the FeCo tip, result in Bi = 17 mT. This larger 
measured value results from the small increase in Bi experienced as one approaches the RF 
microwire source. 

We calculate the Rabi frequencies at Bi = 16 mT for the remaining isotopes to be 680 
kHz for ^H, 165 kHz for ^^Ga, and 117 kHz for "^^As. All isotopes satisfy the adiabaticity 
condition a ^ 1. Due to cubic symmetry, no quadrupolar splitting should be present in 
crystalline GaAs. A small amount of strain due to the mounting process can result in a 
non-zero quadrupolar frequency Qq, though this splitting is likely to be on the order of 
10 kHz for all isotopes. In addition, for nuclear sites near the surface of the particle 
where symmetry is broken, electric field gradients can result in large quadrupolar splittings. 
Nevertheless, the large majority of the nuclear spins detected in our experiments satisfy 
P ^ 1. These two conditions should allow our frequency sweep waveforms to adiabatically 
invert all four nuclear species. 

VII. MODEL AND ESTIMATES 

Modeling the magnitude and shape of the resonance peaks shown in Fig. 4, requires both 
knowledge of the spatial dependence of -Btotai(^) and knowledge of the shape and position of 
the sample. Since -Btotai(0 is strongly inhomogeneous, there is a specific region in space at 
which the magnetic resonance condition is met for each /rf- Only spins near these positions 
are adiabatically inverted and therefore included in the MRFM detection volume. This so- 
called "resonant slice" is a shell-like region in space above the magnetic tip whose thickness 
is determined by the magnetic field gradient and the modulation amphtude fl^p / (27r) of the 
frequency sweeps. We can model this region more exactly using an effective field model for 
adiabatic rapid passage in the manner of Section 4 of the supporting information in Degen 
et al.^ This model shows that the spatial extent of the resonant slice can be described using 
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FIG. 5: Nutation measurement for ^^Ga at T = 1 K. from ^^Ga spins is measured as a function 
of pulse length. A rotating-frame RF magnetic field amplitude Bi of 16 mT is obtained from a 
decaying sinusoidal fit (shown in red) to the Rabi oscillations. Error bars represent the standard 
error of calculated as in Degen et al.^^ 

a simple function: 

^(^-.) ^ _ 7B....ig-2./aF j fo^ (75total(r) " 27r/Rp) < ^Rp 

//(r) = for (75totai(r) - 27r/RF) > ^^rf- (5) 

?7(r*) is normalized to 1 for a nuclear spin positioned exactly in the middle of the resonant 
slice (7-Btotai(0 = '^t^Irf), signifying that this spin is fully flipped by the adiabatic passage 
waveform and contributes its full force to the MRFM signal. A slightly off-resonant spin 
with 1 > ?7(r) > is partially flipped and contributes a fraction of its full force to the MRFM 
signal. Spins outside the resonant slice with ri{r) = contribute no signal. 

In order to calculate the ap, we must therefore determine the intersection of the resonant 
slice with sample for each /^p- In addition, since the gradient varies throughout the resonant 
slice, equal numbers of nuclei at different positions in the slice contribute different forces to 
the final signal. Using (2), (4), and (5) we can then integrate over the volume of the sample 
to find the total MRFM force variance: 
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where S is the sample volume and A is a constant - usually close to 1 - which depends 
on the correlation time of the statistical spin polarization and the measurement detection 
bandwidth. 

We determine -Btotai(^ using a method employed in other recent MRFM experiments. ^'^^ 
First, we measure -Btotai at several different positions above the FeCo tip. The maximum 
value of /rf for which a signal is obtained corresponds to the frequency where the 
resonant slice barely intersects hydrocarbon surface layer closest to the FeCo tip. At this 
frequency /Rp.max, 5totai(n)) = y^RF.max where tq is the position of the hydrocarbon layer 
closest to the FeCo tip. Several such measurements of -Btotai at different tq are then used 
to calibrate a magnetostatic model of the FeCo tip. We infer the shape of the FeCo tip 
from SEM images and we assume a magnetization of 10^ A/m as in previous works. ^'^^ The 
geometrical parameters are fine-tuned in order to produce a field profile Btip{f) which agrees 
with the measured values of -Btotai(^o) = |B + Btip(ro)| foi' our known applied field B. Our 
approximate model then gives us the ability to calculate both -Btotai(0 and ^^^|^(r) at any 
position f. 

Given our approximate knowledge of the shape of the sample from SEM images such as 
Fig. 1, we can only estimate the sample volume S. The GaAs particle is modeled as a 2.4 
/im X 0.6 fim x 0.1 fim rectangular solid with a 200-nm-thick layer of Pt on the end-face. 
The hydrocarbon layer is modeled as a thin film on the surface of this solid. The dimensions 
of this sample are meant to match the cross-sectional size of the particle closest to the FeCo 
tip since this part of the sample contributes nearly all of the observed ap. The back part of 
the sample, with larger cross-sectional area, contributes a vanishingly small ap due to the 
rapid decrease in as a function of distance from the FeCo tip. 

We then use our models for Btotai(0 and S together in a numerical integration of (6) to 
calculate the dependence of ap on /rf- As shown in Fig. 6, the model reproduces the ex- 
perimental data despite the approximate knowledge of the sample shape. Detailed structure 
within the peaks, however, is impossible to reproduce as it is often due to the nanometer- 
scale morphology not included in our idealized geometries. In fact, we can reproduce such 
large variations in the resonance peak shape by altering the details of the sample geometry 
used in our calculation. Prominent structure is particularly evident in resonances measured 
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FIG. 6: MRFM signal of the statistical polarization of (A) ^^Ga and (B) from Fig. 4 along 
with the corresponding calculated signal from the our MRFM model. The sample position used in 
the model was matched to the experimentally set position within an error of 50 nm, likely due to 
experimental position drift during these long time scans. The position was fine-tuned within this 
error range to best match the measured signal. 

with small tip-sample spacings, where the magnetic field gradients are largest and small 
volumes of spins can contribute large force variances. 

A thickness of 2 nm is chosen for the hydrocarbon layer in our model in order to produce 
a resonant approximating our measurements. This thickness is consistent with previous 
measurements of such layers which estimated a thickness of approximately 1 nm.'^'^^ The 
small discrepancy could be due to differences in the surface properties of our sample including 
roughness and affinity to adsorption of hydrocarbons compared to previous samples. 

Despite the approximate nature of our model for ap, we can use it to make an order of 
magnitude estimate of the detection volume in our experiments. Using the parameters of 
each measurement, we can estimate the detection volume Vd as the sample volume intersect- 
ing the resonant slice, i.e. the volume in which (7-Btotai(0 — '^'^Irf) < ^rf- The number of 
spins contained therein is then = naVd- In the case of the peak a'^p from the hydrocarbon 
layer at /rf = 115.5 MHz in Fig. 6, we calculate a = (40 nm)^ and A^^ = 6 x 10^. For this 
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spin ensemble the ratio of SNP to BNP is 0.20. Furthermore, we can estimate the sensitivity 
of this measurement since we know that SNR of our measurement increases with the square 
root of the averaging time. We calculated the SNR at each /rf by dividing the measured 
(jp by the standard error of this measurement calculated as in Degen et al.^^ This error 
takes into account both the noise due to fluctuations in the cantilever motion, i.e. thermal 
noise and non-contact friction, and the noise due to the statistically polarized spin ensemble 
itself. Given the SNR of 14.4 achieved after 300 s of averaging, we estimate a measurement 
sensitivity equivalent to 7 x 10® spins/v^Hz. In general, the sensitivity of these measure- 
ments is limited by the mechanical fluctuations of the cantilever due to thermal noise and 
non-contact friction. 

We can make similar calculations for the quadrupolar nuclei. The peak values of ap shown 
in Fig. 4, however, do not represent the maximum attainable signal for each isotope. Due 
to the long averaging times required for these isotopes, position scans used to optimize the 
signal amplitude were not performed before these measurements. Approximate measurement 
positions were estimated based on the experiments resulting in smaller than optimal a p. 
For ^^Ga, however, an x and y position scan was performed in order to find the optimal 
ap = 25 aN^ at /rf = 34.95 MHz. From this scan it was found that changes in position less 
of only 50 nm resulted in signal loss of over a factor of 2, emphasizing the importance of 
optimal alignment. This signal is in reasonable agreement with the maximum signal ap = 20 
aN^ calculated for the same parameters in our model. Using our model we can calculate 
Vd = (260 nm)^ and = 2 x 10^ for the ''^Ga signal plotted in Figs. 4 and 6, which is 
slightly shifted from the optimal position. In this case, we find SNP/BNP = 0.14. The 
corresponding sensitivity is estimated from the SNR of 15.4 after 600 s of averaging to be 
2 X 10*^ "^^Ga spins/v^Hz. Similar calculations are not carried out for ®^Ga and ''^As, though 
sensitivity for these isotopes should be of the same order after a scaling factor equivalent to 
appropriate MRFM receptivity. 

As discussed in Section V, the large difference in the sensitivity between and the 
quadrupolar nuclei is mostly due to the 200 nm Pt layer which forces the Ga and As moments 
into a region of far smaller magnetic field gradient than at the hydrocarbon layer. Future 
experiments should be designed such that this Pt layer, which is an artifact of the FIB 
mounting process, is not present. Without this intermediate layer, far better sensitivities 
should be achieved for the quadrupolar nuclei. Table III shows predicted sensitivities for 
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iH 


7iGa 




(hydrocarbon layer) 


(GaAs) 


Sample-tip distance (nm) 


100 


300 


Maximal at sample (T/m) 

ox f \ 1 J 


5 X 10^ 


8 X 10^ 


Maximal Btip at sample (T) 


0.10 


0.03 


/rf (MHz) 


115.5 


34.95 


ORF/(27r) (kHz) 


400 


100 


Nd 


6 X 10^ 


2 X 10^ 


Vd 


(40 nm)^ 


(260 nm)3 


SNP/BNP 


0.29 


0.14 


Averaging time (s) 


300 


600 


Sensitivity (spins/\/Hz) 


7 X 10^ 


3 X 10® 



TABLE H: Detection and sensitivity estimates for the H and Ga resonances plotted in Fig. 4 at 
B = 2.65T and T = IK. 

^^Ga, '^^Ga, and ^^Ga for a 100 nm spacing between the sample and the FeCo tip - without 
any intermediate layer. All other parameters are identical to those of the actual experiments. 
These extrapolations are based on positioning the Ga and As in the same position as the 
nuclei in our experiment. We make the assumption that the noise would be the same as 
that measured in the experiment. 

VIII. CONCLUSION 

The results presented here, demonstrate our ability to detect nanometer-scale volumes 
of Ga and As nuclei using MRFM. Given the spin sensitivity extrapolated from our data 
and our model, the detection of III-V nanostructures such as nanowires or sub-surface self- 
assembled InAs QDs should be possible. Self-assembled InAs QDs, for example, contain 
10^-10'' nuclear spins, lie as close as 50 nm from the wafer surface, and could be attached 
to a cantilever using the FIB technique demonstrated here. Further improvements to the 
force sensitivity - most importantly for reducing measurement times - will be required in 
order to realize MRI in III-V materials with better than 100 nm resolution. The potential 
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69Ga 


7iGa 


75 As 




(GaAs) 


(GaAs) 


(GaAs) 


Sample-tip distance (nm) 


100 


100 


100 


Maximal at sample (T/m) 

CfX f \ 1 J 


5 X 10^ 


5 X 10^ 


5 X 105 


Maximal Btip at sample (T) 


0.10 


0.10 


0.10 


/rf (MHz) 


27.9 


35.1 


19.8 


ORF/(27r) (kHz) 


100 


100 


50 


Nd 


1 X 10^ 


9 X 10^ 


2 X 10^ 


Vd 


(210 nm)3 


(210 nm)3 


(200 nmf 


SNP/BNP 


0.20 


0.19 


0.23 


Calculated a], (aN^) 


250 


210 


150 


Sensitivity (spins/\/Hz) 


1 X 10^ 


9 X 10^ 


3 X 10^ 



TABLE HI: Extrapolated detection and sensitivity estimates for the quadrupolar nuclei based on 
parameters achieved for ^H at 5 = 2.65T and T = IK in Table II. 

for sub-surface, isotopically selective imaging on the nanometer-scale in III-V materials is 
a particularly exciting prospect since conventional methods such as SEM and TEM lack 
isotopic contrast. 
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where A is a measurement bandwidth much smaller than the cantilever bandwidth. 
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